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ABSTRACT 

q ' We present an analysis of brown dwarf model spectra in the mid-infrared spec- 
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tral region (5 - 20 /xm), in anticipation of data obtained with the Space Infrared 
Telescope Facility. The mid-infrared spectra of brown dwarfs are in several ways 
simpler than those in the near-infrared and yet provide powerful diagnostics of 
brown dwarf atmospheric physics and chemistry, especially when combined with 
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CSj ; ground-based data. We discuss the possibility of detection of new molecular 

species and of the silicate cloud, predict strong observational diagnostics for non- 
equilibrium chemistry between CO and CH 4 , and N 2 and NH 3 , and speculate on 
^ ■ the possibility of discovering brown dwarf stratospheres. 
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Subject headings: stars: low mass, brown dwarfs — radiative transfer — molec- 
ular processes — infrared: stars 



Introduction 



The near-infrared and optical spectra of brown dwarfs have received intense observa- 
tional and theoretical scrutiny (see the recent reviews by Basri 2000 and Burrows et al. 
2001) but their mid-infrared spectra have been relatively neglected both in observation and 
theory. There are only two announced photometric detections of a brown dwarf beyond 5 /im 
(Matthews et al. 1996; Creech-Eakman et al. 2003) and neither strongly constrains models. 
Although synthetic mid-IR spectra have been published (e.g. Marley et al. 1996; Tsuji, 
Ohnaka & Aoki 1999; Allard et al. 2001; Burrows, Sudarsky & Lunine 2003) there has yet 
been no detailed discussion of the expected spectral diagnostics in L and T dwarfs. With 
the Space Infrared Telescope Facility (SIRTF) now in orbit, the dramatic improvement in 
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our ability to observe brown dwarfs in the mid-infrared will give new insights in the astro- 
physics of the complex atmospheres of these cool, dim denizens of our neighborhood. In 
this Letter, we present mid-IR model spectra of brown dwarfs and discuss their properties 
in terms of effective temperature, important molecular absorbers, the role of silicate clouds, 
and anticipate potential discoveries with SIRTF. 

The instruments onboard SIRTF cover the wavelength range from 3 to 180 /im. Two 
of these will target brown dwarfs: the Infrared Array Camera (IRAC) and the Infrared 
Spectrograph (IRS). Brown dwarfs will be imaged by IRAC in four bandpasses centered 
at 3.6 ^m, 4.5 fim, 5.8 jj,m, and 8.0 /um, respectively (Fig. I) 4 . Since brown dwarfs become 
very dim and their spectra contain little information beyond 20 yum, the most useful IRS 
observations will be in the "Short wavelength, Low resolution" (SL) and the less sensitive 
"Short wavelength, High resolution" (SH) modes 5 . The SL and SH modes cover the 5.3 - 
14.2 /Ltm and the 10.0 - 19.5 (im spectral bands, respectively. 



2. Atmosphere models and spectra 

To model the atmospheres and spectra of the L- and T-dwarfs we employ the radiative- 
convective equilibrium atmosphere model of Marley et al. (1996; further described in Burrows 
et al. 1997 and Marley et al. 2002) which includes the precipitating cloud model of Ackerman 
& Marley (2001). High resolution spectra are computed from the initial temperature profile 
and cloud structure (Saumon et al. 2000; Geballe et al. 2001). The chemistry is computed 
in the framework of the cloud condensation model (Lodders 1999a; Lodders & Fegley 2002; 
Lodders 2002; Lewis 1969). The opacity includes the molecular lines of H2O, CH4, CO, 
NH 3 , H 2 S, PH 3 , TiO, VO, CrH, FeH, C0 2 , HCN, C 2 H 2 , C 2 H 4 , C 2 H 6 complemented with 
the atomic lines of the alkali metals (Li, Na, K, Rb and Cs) and continuum opacity sources 
from H 2 CIA, H 2 , H and He Rayleigh scattering, H~ bf and ff, H 2 ff, He~ ff, and H 2 bf and 
ff. While the chemical equilibrium is computed with a large number of condensates, only 
Fe, MgSi0 3 , A1 2 3 , H 2 0, and NH 3 are considered in the cloud model. Condensed Mg 2 Si0 4 
is accounted for with MgSi0 3 and the remaining condensates are not appreciable sources of 
opacity. 

We consider spectra from T eff = 600 to 2400 K which covers the full range of L and 
T dwarfs. We limit the discussion to models with solar metallicity. The sedimentation 
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parameter of the cloud model is fixed at / sc d = 3, which gives a good representation of far- 
red and near-IR photometry of L dwarfs (Marley et al. 2002; Burgasser et al. 2002) and of 
the ammonia cloud deck of Jupiter as well (Ackerman & Marley 2001). For simplicity, we do 
not take into account the possibility of cloud disruption near the L/T boundary (Burgasser 
et al. 2002). 

3. Mid-infrared spectra of brown dwarfs 

The combination of IRAC photometry and IRS spectroscopy with ground-based op- 
tical and near-IR data will give the complete spectral energy distributions (SED) of many 
brown dwarfs. Empirical bolometric corrections immediately follow as well as the bolometric 
luminosity for objects with known parallaxes. An independent determination of T e fj (e.g. 
by fitting spectra with models) will give the gravity, radius, and mass of individual brown 
dwarfs. 

It is well known that in the far-red and near-IR, the brightness temperature Tb r of brown 
dwarfs varies strongly with wavelength due to the high contrast between opacity windows and 
molecular absorption bands. This is still true in the mid-IR up to about 11 [im, thereafter 
Tb r gradually decreases to stabilize at ~75 % of T e g beyond 18 //m. Only near 6 ^m and 
near 10 /xm does T br become as large as T eS . For 4 < logg(cgs) < 5.5, the pressure at 
the mid-IR photosphere stays between 0.1 and 3 bar, depending mostly on wavelength and 
surface gravity and least on T c s. Basically, the mid-IR spectra of brown dwarfs are formed 
at T ^ T c s and pressures of about 1 bar. 

3.1. Molecular species 

Figure 1 identifies the mid-IR molecular absorbers in a sequence of spectra from 600 to 
2400 K. By far, the most important absorber is H 2 throughout the entire spectral range, 
except for two strong molecular bands due to CH 4 and NH 3 . Beyond 15 ^m, all features are 
due to H 2 0, with the exception of a few weak NH 3 bands between 32 and 46 /zm for models 
with T c ff ^ 800 K. In the high-T e ff spectra, the fundamental band of CO at 4.8 /im and a weak 
TiO band at 10 fim are the only features not originating from H2O. The step in the spectrum 
~ 6.5 /im is a H 2 feature. The TiO band disappears at 2200 K but the CO band persists 
down to 900 K in these cloudy models. As T e s decreases, new molecular bands appear and 
steadily increase in strength. The 10.5 fxm band of NH 3 appears at 1150 K. This band is very 
broad and dominates the spectrum from 8.5 to 16 at T eS < 800 K. Weaker bands of NH 3 
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appear at 900 K (5.5 - 7 /xm) and 800 K (3.9 - 4.5 /xm). The former is superimposed on a 
strong H 2 band, however. Methane shows only one band in the mid-IR, centered at 7.8 /im. 
It appears at 1600 K and becomes very strong at low T e g to dominate the spectrum between 
7 and 9.2 /im. Broad, featureless CIA opacity appears faintly in the 9 to 13 /im region in 
cool models below 1000 K, but it shows no gravity sensitivity and is probably undetectable. 

Our current phosphorus chemistry uses updates for several P-gases (PH, PH3, PN, PS; 
Lodders 1999b, 2004) and the P4O6 gas thermodynamic properties from Gurvich, Veyts & 
Alcock (1989) instead of those from Chase (1998). Our previous models used P4O6 data from 
Chase (1998) and calculated spectra did not show the strongest PH 3 band centered at 4.3 /im, 
just to the blue of the fundamental CO band. However, the Gurvich et al. data make Pzt0 6 
significantly less stable which results in a larger abundance of PH 3 and the band appears 
below HOOK with a characteristic narrow absorption spike. This feature should provide a 
good observational test of the phosphorus chemistry in brown dwarfs, but its signature in 
IRAC band 2 photometry will have to be distinguished from possibly enhanced CO (Section 
5). 

The most abundant sulfur-bearing compound is H 2 S, a molecule that has yet to be 
detected in brown dwarfs. Mid-IR observations will not improve its prospects of being 
noticed since its opacity is consistently 3-4 orders of magnitude below that of H 2 0. 

We predict that CO and CH 4 will be visible simultaneously in the mid-IR for T cff ^ 
1600 K, i.e. from mid-L to mid-T spectral types. The NH 3 band at 10.5 /xm will be easy 
to detect and is a characteristic of all T dwarfs (see Section 5, however). Observations of 
Gl 229B should confirm the weak detection of this molecule in the near-IR (Saumon et al. 
2000). 

Figure 1 shows that the mid-IR spectra of brown dwarfs are characterized by molecular 
bands with strong T e g dependences below 1600 K. Fitting the IRS spectra alone with models 
should allow the determination of T e g of brown dwarfs of mid-L and later spectral types. 

3.2. Silicate cloud 

The effects of silicate, iron, and corundum clouds on the near-IR spectra of brown dwarfs 
have been discussed extensively (Ackerman & Marley 2001; Allard et al. 2001; Burgasser 
et al. 2002; Cooper et al. 2003; Marley et al. 2002; Tsuji 2002; Tsuji & Nakajima 2003). 
Because the mid-IR photosphere is characterized by a fairly constant brightness temperature 
and a single dominant absorber (H 2 0), the spectral signature of clouds stands out more 
clearly than in the near-IR. Most interesting is the possibility of detecting the silicate cloud 
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opacity feature ~ 10 jum. Figure 2 shows spectra computed with and without the cloud 
opacity using the same cloudy atmospheric structure. We find that the silicate feature (a 
rise in the opacity between 9 and 10 /xm to a new plateau) is masked by a similar behavior 
in the H2O and CH4 opacities. Nevertheless, the cloud opacity is large enough to appear 
in two windows where the gas is less opaque in models with 1200 ^ T c s ^ 2000 K: from 
~ 4 to 6.5 /Ltm (affecting IRAC bands 1 through 3) and from 9 to 13 fim. The effect is 
strongest at T eS ~ 1600 K, which corresponds to the largest cloud optical depth above the 
10 fxm photosphere. Because the silicate absorption flattens the spectrum noticeably in the 
9 to 12 /Ltm range, its identification should be fairly straightforward in IRS spectra of mid-L 
dwarfs. For models with a higher condensate sedimentation efficiency (/ se d = 5), the thinner 
cloud layer remains optically thin above the photosphere and is invisible at 10 jj,m. This 
feature can thus serve as a diagnostic of the vertical structure of the cloud model as its 
strength depends strongly on the cloud's vertical extent and on the particle size. 



3.3. IRAC colors of brown dwarfs 

The IRAC will image brown dwarfs and provide photometry in 4 broad bandpasses (Fig. 
1), labeled Band 1 through 4. Bands 1 and 2 cover wavelengths not accessible to the IRS. 

As can be seen in Fig. 1, Band 2 covers the bands of CO and PH 3 but is nevertheless the 
least affected by strong molecular bands over the full range of T eff of interest. Bands 1 and 
4 are increasingly affected by CH 4 absorption below T c s ~ 1500 K. Strong H 2 absorption 
appears in Band 3 below T cS ~ 1200 K. 

A study of color-color diagrams shows that IRAC colors will be most useful as diagnostics 
of brown dwarf physical properties for T eff ^ 1400 K. At higher T eS , the colors change little 
with T cff or the T eS and gravity dependences become hopelessly tangled. For brown dwarfs 
with known distances, color-magnitude diagrams provide a good discriminant of both T cS 
and gravity. 



4. Trace molecular species 

The IR spectra of giant planets reveal the presence of trace carbon species such as C0 2 , 
C2H2, and HCN in emission. These arise primarily from the photochemistry of CH4 (also 
seen in emission) in the stratosphere, i.e. above a temperature inversion in the atmosphere 
(Moses 2000). No published brown dwarf model shows a temperature inversion, mainly 
because the present knowledge of the physics of brown dwarf atmospheres is too primitive to 



5 



model a credible stratosphere. We nevertheless anticipate that brown dwarfs also may have 
stratospheres (Yelle 2000) and observations in the mid-IR are a powerful way to discover 
them by detecting emission features from trace molecules such as CO2, C2H2, C2H4, C2H6, 
and HCN that should arise from photochemistry driven by the stratospheric and possibly 
the weak interstellar UV flux. In L and T dwarfs, H 2 remains in the gas phase and can 
participate in the photochemistry, leading to the formation of species such as HCO and 
CH 2 (Friedson, Wilson & Moses 2003). 

In the absence of a stratosphere, these molecules would appear in absorption. The 
equilibrium abundances, i.e. in the absence of photochemistry, are all quite low. We have 
calculated the minimum enhancement factor (e) of the abundance of several trace species 
required for detection. We define e as the ratio of the column density of the trace species 
above the photosphere that is needed for detection to the equilibrium column density above 
the photosphere. For each species in Table 1, we give the wavelength of the strongest opacity 
feature, the T e g- that is most favorable for detection (i.e. giving the lowest enhancement 
factor), and the enhancement factor. Clearly, none of these species is likely to be detected in 
absorption, except for CO2, which would be favored by higher metallicity and lower gravity. 
For instance, the equilibrium abundance of CO2 should be marginally detectable at 15 ^m 
in a model with T c g = 1200 K and log g = 4. 

Objects with masses below 13 Mj do not burn deuterium and deuterated molecular 
species are expected in their atmospheres. With its distinctive 4.55 ^m band and a relatively 
large abundance, CH3D is the most easily detected deuterated species expected in brown 
dwarfs. We find that for CH3D/CH4 = 2 x 10~ 5 , an enhancement factor of 250 would be 
necessary to detect it in T dwarfs above T eS = 600 K (Table 1). Detection becomes easier at 
lower T cS but is very unlikely among the known objects. 

5. Non-equilibrium chemistry 

Vertical transport in brown dwarf atmospheres can lead to non-equilibrium abundances 
in slowly reacting species such as CO and N 2 . The net effect is to enhance the upper 
atmosphere abundances of CO and N 2 to the detriment of CH 4 , H 2 and NH 3 (Fegley & 
Lodders 1996; Griffith & Yelle 1999; Lodders & Fegley 2002). The detection of greatly 
enhanced CO (Noll, et al. 1997) and of depleted NH 3 (Saumon et al. 2000) in Gl 229B, as 
well as the M' flux deficiency observed in several T dwarfs (Golimowski et al. 2004) suggest 
that non-equilibrium chemistry is a general phenomenon in brown dwarfs. Nevertheless, it 
remains mostly unexplored observationally. 
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We have modeled the non-equilibrium spectra of brown dwarfs following the method 
described in Saumon et al (2003). We find that the 7.8 fim band of CH4 appears at T e g ~ 
1400 K and 1200 K for mixing coefficients of K zz = 10 2 and K zz = 10 4 cm 2 /s, respectively, 
compared to 1600 K in the case of equilibrium chemistry. The 10.5 fj,m band of NH3 appears 
at T c g ^ 1000 K for K zz ^ 10 2 cm 2 /s, instead of 1150 K. Because of the the relatively flat 
equilibrium abundance profile of NH 3 , this result is nearly independent of K zz . While NH 3 
can be depleted by more than one order of magnitude (for T cS = 800 K), the 10-11 /xm 
features remain strong (Fig. 3). 

6. Conclusions 

Spectroscopy with the IRS instrument on SIRTF between 5 and 20 /j,m will complete 
the sampling of the spectral energy distribution of brown dwarfs. From our analysis of our 
mid-IR synthetic spectra, we find that the new data will: 1) clearly reveal the presence of 
NH 3 in T dwarfs, despite the strong depletion due to vertical transport in the atmosphere, 
2) likely detect the silicate cloud ~ 10 /im in mid-L dwarfs and put strong constraints on the 
vertical structure of the cloud, 3) lead to a much more complete picture of non-equilibrium 
chemistry among CO, CH4, H2O, N2 and NH3, 4) show CO2 in low-gravity, high metallicity 
targets near the L/T transition, 5) not reveal interesting species such as H 2 S, CH 3 D, and 
the H2 CIA opacity, and 6) possibly discover brown dwarf stratospheres through emission 
from trace species such as CO2, HCN, HCO, C2H2, C2H4, C2H 6 and CH 2 0. Furthermore, 
we predict that: 7) a strong band of PH 3 at 4.3 /jm falls in the IRAC band 2, 8) IRAC 
photometry will be most useful for objects below 1400 K, and 9) photospheric spectra of 
brown dwarfs are of very little interest beyond 20 fxm. The mid-IR spectral window is rich 
in diagnostics of brown dwarf atmospheres. Our understanding of these cool neighbors will 
grow dramatically with our first mid-IR observations with SIRTF. 

We thank Richard Freedman for providing the molecular opacities used in this cal- 
culation and the SIRTF/IRS instrument team for informative discussions. This work was 
supported in part by NASA grants NAG2-6007 and NAG5-8919 and NSF grant AST 00- 
86288 (MSM) and by NSF grant AST 00-86487 (KL). Part of this work was supported by 
the United States Department of Energy under contract W-7405-ENG-36 
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Fig. 1. — Sequence of mid-IR spectra of brown dwarfs with T eff ranging from 600 to 2400 K 
(from bottom to top) in steps of 200 K. All models have solar metallicity, log g — 5, / se d = 3 
and are shown at a spectral resolution of 200. Spectra are shifted vertically for clarity. 
Molecular bands that appear at high T e g are identified at the top, those that appear at low 
T e fj are shown at the bottom. Outside of those bands, H 2 is responsible for all features. 
Weak bands are indicated with dashed lines. The bandpasses of the four IRAC filters are 
shown by thick solid lines at the top. 

Fig. 2. — Effect of silicate cloud opacity on the mid-IR spectra of brown dwarfs. The spectra 
shown have T c g = 1000 to 2400 K in steps of 200 K. The dotted lines show spectra computed 
from the same cloudy (T, P) structures and chemistry, but without the cloud opacity. IRAC 
bands 3 and 4 are shown by thick solid lines at the top. 

Fig. 3. — Effect of non-equilibrium chemistry on the mid-IR NH 3 band. The spectra shown 
have T e g = 800 to 1600 K in steps of 200 K. The spectra shown are computed computed from 
the same cloudy (T, P) structures but with equilibrium (solid) and non-equilibrium (dotted) 
chemistry. See Fig. 2 for details. 
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Table 1. Enhancement factors for trace species. 



Species 


A(/xm) 


T eff (K) 


log(e) 


C0 2 


15.0 


1400 


0.7 


HCN 


14.0 


1400 


2.7 


C2H2 


13.7 


1400 


5.8 


C2H4 


10.5 


1200 


4.9 


C2H 6 


~ 12.2 


1000 


3.5 


CH 3 D 


4.55 


< 600 


< 2.4 b 



a For log g = 5 and solar metallicity. 
b Cloudless atmosphere model. 
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